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For oscillation phenomena which may also have an unstable course, 
it is customary to represent the regions where stability or instability 
prevails in a plane as functions of two parameters x and y. 

In order to determine whether stability exists at any point of the 
plane represented (thus, whether a disturbance of the oscillation phe- 
nomenon considered is damped in its course) it must be investigated 
whether all roots X of the frequency (characteristic) equation 


X n + apX n “ 1 + . . ..+ 8^ = 0 


have a negative real part at this arbitrary point. The mathematical 
condition for this is known to be that the n Hurwitz determinants Dp 
to Dq which are formed from the coefficients ap to a n and are 

functions of the two parameters x and y mentioned before are all 
positive for this point. Resulting from this criterion and completely 
equivalent to it is the fact that the coefficients a^ to a^ and only 

a few certain Hurwitz determinants must be positive. 

If, conversely, all roots X of the frequency equation have a nega- 
tive real part, all values Dp to 1^ are positive. If one now visu- 
alizes the point considered before (for which we assume stability to have 
been established) as traveling in the representation plane of figure 1, 
there vary with Its parameters -x and 'y also the coefficients ap 

to a n , the real parts of the roots X, and finally the n Hurwitz- 
determinants . If one arrives at a point where for instance a real root 


*"Zur Par ste Hung des Stabilitatsgebietes bei Schwingungsaufgaben 
mit Hilfe der Hurwitz- Determinanten. " Schweizer Archiv, March 1930 , 
pp. 93 - 96 . 
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disappears , a n also disappears since. |a n | is the product of the values 

of all roots; if one reaches, in contrast, a point where -the real part 
of a complex root becomes zero, it can be shown that then Dq _1 becomes 
zero. For every case, however, the product anD n _i disappears which is 
nothing else but the Hurwltz determinant of the nth degree 


a n D n-l = I) n 


Thus the important theorem, the proof for which will be presented 
later, is valid: 

The limits of the (usually unique) stable region lie at 
®n = a n®n-l = 

«- 

For graphical representation, it is therefore completely sufficient 
to plot D n = 0 or more simply a n = 0 and D n _p =0 as separate 

limiting curves of a region for which stability is known to - prevail at 
an arbitrary point, as illustrated by figure 1. If the limit D n _i = 0 

is exceeded, the course of the. oscillation process is "dynamically" 
unstable because a damping becomes negative; beyond the limit a n = 0 , 
one usually calls the oscillation process "statically" unstable. 

In particular, the following is valid for oscillation phenomena 
which lead to frequency equations of the 4th degree (reference l). 

The Hurwitz determinant of the 4th degree formed from the coeffi- 
cients A to E of the frequency equation 


AX 1 *- + BX 3 + CX 2 + DX + E = 0 

reads 


= ED 3 = E(BCD - AD 2 - B%) = ER 


with the expression in parentheses known as Routh's discriminant R; the 
latter is therefore nothing else but the Hurwitz determinant of the 
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3rd degree. Since the coefficient A is usually +1, it is generally 
valid as dynamic stability condition that the expression (BC - D)D - B% 
turns out positive. The static stability is then guaranteed in addition 
by E > 0 so that the graphic representation of the region of figure 2 
results . 

Therewith, every requirement has been met; for it is impossible 
that within the region denoted as stable a curve C = 0 or D = 0 
could run and perhaps still further reduce this region. 

It is therefore completely superfluous to investigate further what 
sign the other coefficients of the frequency equation have once the 
boundaries for a region E = 0 and R = 0 have been fixed; however, 
and this is important, it must be known that, at an arbitrary poipt of 
this region, stability actually prevails. For this it is sufficient 
to check for a point conveniently situated (for instance on an axis) 
the signs of the remaining coefficients A to D of the frequency 
equation which is equivalent with the fact “that there all Hurwitz deter- 
minants Dp to Dp turn out positive. That this is important can, for 
instance, be recognized from the fact that, in the region indicated on 
the left in the plane of representation of figure 2, instability may 
prevail in spite o.f E > 0 and R > 0. — 

The advantage of these recognized facts given above does not result 
In much simplification for frequency equations of the 4th, 5th, or even 
6th degree; however, this advantage may save a great deal of unnecessary 
calculations. The expressions for the Hurwitz determinants become more 
and more cumbersome for higher degrees so that It will probably be a 
very welcome facilitation for the outlining of the stability region if 
the mere representation of the two curves a n = 0 and Dn_i = 0 will 
be sufficient. 

There now follows the proof of the theorem given before that the 
limits of the (usually unique) stable region lie at D n = a^Dj^p = 0. 

1. For the present consideration, the existence of at least one 
stable region is presupposed. Not even one need necessarily always 
exist; on the other hand, s^Peral stability regions, separated from 
one another, may occur. One can easily make sure of this by considering 
a frequency equation of the 2nd degree which pertains to an ordinary 
oscillation 
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The coefficients a^ and a 2 are dependent on two parameters x 
and y, corresponding to the two coordinate axes of the plane of repre- 
sentation. For simplicity's sake, we shall assume that a 2 everywhere 
has positive values only; then the sign of a-]_ alone is the criterion 
for the stability, since aj_ is a measure of the damping of the oscil- 
lation considered and this oscillation is Btable only if aj_ is posi- 
tive. The limit of stability now depends entirely on the sign of the 
analytical (or empirical) function of the coefficient a-j_ of the two 

parameters x and y, or, in other words, what form the curves a-j_ = 0 
have in the plane of representation which separate the positive from the 
negative values. In this manner, one can easily understand that several 
stability regions may exist but that just as well not even one need exist 
anywhere In the entire plane of representation. 

For frequency equations of higher degree, the probability decreases 
that simultaneously In several regions all stability conditions found by 
Hurwitz will be satisfied. Usually, one will deal with only one single 
stability region which will then be considered more thoroughly. 

2. Hurwitz' criterion (reference 2) signifies that- a prescribed 
equation of the nth degree with real coefficients 


a 0*» + a^X 11 "^ + . . . + a n = O 


(&0 >0 ) 


possesses roots exclusively with negative real parts only when the values 
of the n-determinants 



a l 

a 3 

®5 ‘ 

• •" a 2k-l 


Sq 

a 2 

a-ii 

• • s 2k-2 

Dj, = 

0 

a l 

a 3 * 

* * a 2k-3 


• 

• 

• • 

- * a k 


(k = 1,2, ... ,n) 

are all positive. Therein, one generally has to put a x = 0 when the 
subscript x is negative or larger than n. 
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It shall now he proved that the Hurwitz determinant of the (n - l)th 
degree always disappears when the real part of a complex pair of 

roots is zero. 

Following the way Hurwitz pursued in deriving his criterion,, one 
finds the remarkable presupposition that no purely imaginary roots must 
exist. This, however, is precisely the condition whose influence on a 
certain Hurwitz -determinant is of special interest. However, for that 
reason, one need not abandon this presupposition; it is sufficient to 
interpret I> n _i simply as a prescribed calculation rule which one applies 

without consideration of its derivation, to the coefficients of an equa- 
tion of the nth degree. 

It will now be expedient to consider (corresponding to the assump-“~ 
tion that the real part of a complex pair of roots is to be zero) the 
equation of the nth degree as split into two factors: one Is a quad- 

ratic expression with the purely imaginary pair of roots and the other 
Is an expression of the (n - 2)th degree, thus — ■ — 


(l 2 + a)(a 0 A . n “ 2 + a-L *. 11 " 3 + 


+ a n-2 



0 


therein one may put the coefficient Sq = 1, without impairing the gen- 
erality; the designation &q will, however, be retained. For all 
remaining coefficients aj_ to a n _ 2 , no additional restriction will be 

prescribed other than that they are to be real and not all zero (which 
would be trivial for this statement of the problem) . Multiplication of 
the two factors yields 


SqXP + &2^ P' ^ + ^ag + aao)l n 2 + (^3 3 + • • . + 

Aq A 1 A 2 A 3 



The combined coefficients A one now visualizes as substituted in 
the Hurwitz determinant D n t, the general term of which has in the 
i 
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sth column and zth line the subscript x = 2s - z. Thus, one has there 
&2s-z OT > expressed by the coefficients a, the term 


a 2s-z + aa 2s-z-2 


The value of the determinant D n _2_ is now extended by multiplying 
all terms of the first line by Sq and all terms of the second line 
by Furthermore, one makes use of the. theorem that a determinant 

does not change its value if one adds to the elements of one line the 
elements of another line multiplied by an arbitrary number. This is done 
by adding to the 'terms of the first line , which already have been multi- 
plied by Sq, the terms of the third, fifth, . . . line_multiplied by 

a 2 } a V * * •* ttius , one obtains as the general term of- the sth column 
in the first line 



z-1, 3,5, 


« • » 



aa 2s— z— 2 ) 


Likewise, one adds to the terms of the second line already multiplied 
by a-j_ the terms of the fourth, sixth, . . . line multiplied by a^, 

a^, . . . ; one then finds as the general term of the sth column in the 
second line 



z-2,4, 6, 



aa 2s-z-2) 


In order to calculate these sums, it is sufficient to note that one 
must generally put a x = 0 when the subscript x is negative; there 

then results for a term of the first line 
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H 

z— 3 , 5 , . . . 


a z-l( a 2s-z + aa 2s-z-2) “ a 0 a 2s-l + a £ a 2s-3 + 


+ 


a 2s-4 a 3 + a 2s-2 a l + a ( a O a 2s-3 + a 2 a 2s-5 + • • • + a 2s-l*- a l) 


and for a term of the second line 



z=2 ,4,6, 


a z-l( a 2s-z + aa 2s-z-2) " a l a 2s-2 + a 3 a 2s-4 + 



* 


a 2s-3®2 + a 2s-L a 0 + a ( a l a 2s-l<- + s 3 a 2 s -6 + 


+ a 2s-3 a o) 


One can see that both sums are of equal magnitude which., however, 
does not signify anything else hut that the corresponding elements of 
the first two lines are equal and that therefore the Hurwitz-determinant 
D n _-]_ extended by identically disappears: 


ViVi * 0 


From the derivation follows that one could have extended initia lly , 
instead of the first and second line, two arbitrary other odd and even 
lines by a s _-]_ and then have proceeded further in the same manner; one 

then would have obtained quite generally 


a a ,,D , = 0 

even odd n-1 


% 
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Under the obvious assumption that there always exists such a pair 
of values of the coefficients 8 q to a n _2> the Product of which does 

not disappear, thus 


= 0, Q.E.D# 


is valid. 

3. From the derivation of this proof, one may further conclude that 
part of the curve D n _j_ = 0 may run entirely in the unstable region; 

thus it does not appear there as the required stability limit. Since 
such a possibility had been pointed out before, the motivation of this 
noteworthy indication should be mentioned. 

In splitting the prescribed equation of the nth degree into two 
factors, no more detailed data on the coefficients ao to a n _2 of 

the expression of the (n - 2)th degree have been given and accordingly 
none regarding the roots of the equation 



A *- 2 


+ 


,n-l 




• + a n-2 



0 


either. 

These roots may therefore have negative as well as_positive real 
parts in the neighborhood of points of the plane of representation for 
which Dn-1 disappears; in other words, they may signify stability and 
also instability. If there exists at least one root ofhihe former 
expression of the (n - 2)th degree with a positive resl_part, this fact 
signifies that the curve D n _-j_ = 0 lies in an unstable region. 

Figure 3 is to indicate how, for instance, three roots may vary 
their position in the complex number plane If rtheir point of reference 
in the plane of representation shifts beyond the curve P n _ 2 _ = 0. Thus 

D n _p = 0 is not necessarily always the boundary between' two regions 
where stability and instability prevail. 

Neither does the sign of the coefficient ..a in the_quadratic expres- 
sion + a) play a role in the derivation of ..the proof. If a is 
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positive, a purely imaginary pair of roots satisfies the equation 
(A,2 + a ) =0; this case has been considered in particular. However, 
a could just as well be negative without altering the result of the 
theorem. This signifies, however, that in the presence of two real 
opposite -equal roots of the equation of the nth degree, the Hurwitz- 
determinant of the (n - l)th degree also disappears. 

As figure k shows, here also an unstable region lies to both sides 
of the curve D n _i =0 so that again it does not appear as the sta- 
bility limit. 

Thus, these examples show that it is indispensable to make sure 
whether really all stability conditions on one side of the curve 
D n i = 0 have been fulfilled; only then that curve forms together with 

the curve a n = 0 the boundary against the unstable region. 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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Figure 3.- Change of position of the roots in the complex number plane. 
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Figure 4.- Unstable region for. the curve D n _^ =0. 
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